A single-ion nonlinear mechanical oscillator 
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We study the steady state motion of a single trapped ion oscillator driven to the nonlinear regime. 
Damping is achieved via Doppler laser-cooling. The ion motion is found to be well described by 
the Duffing oscillator model with an additional nonlinear damping term. We demonstrate a unique 
ability of tuning both the linear as well as the nonlinear damping coefficients by controlling the 
cooling laser parameters. Our observations open a way for the investigation of nonlinear dynamics 
on the quantum-to-classical interface as well as mechanical noise squeezing in laser-cooling dynamics. 
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Nonlinear dynamics prevails in many dynamical sys- 
tems in nature, introducing a rich behavior such as criti- 
cality, bifurcations and chaos. Nonlinear dynamics on the 
microscopic scale is especially interesting as it can shed 
light on the quantum-to-classical transition as well as 
provide a mean to suppress thermal and quantum noise. 

All mechanical oscillators will show nonlinearity when 
driven far enough from equihbrium. The simplest such 
nonlinear oscillator is the Duffing oscillator which in- 
cludes a cubic term in the restoring force [T]. Recently, 
such Duffing nonlinear dynamics has been extensively 
studied with nano-clcctromechanical beam resonators. 
The basins of attraction of a nano beam oscillator were 
mapped H]. Noise squeezing and stochastic resonances 
were observed close to the Duffing instability [31 U • Noise 
squeezing was predicted to enable mass and force detec- 
tion with precision below the standard thermal limit 0- 
[7] and possibly below the standard quantum limit when 
operating close to the oscillator ground state [8]. 

The mechanical motion of trapped ions is highly con- 
trollable and can be efficiently laser-cooled to the quan- 
tum ground state |9]. High fidelity production of Fock, 
squeezed, and Schrodinger-cat states was demonstrated 
with a single trapped-ion [TUl [TT| . At the temperature 
range obtained with laser-cooling techniques, quadruple 
RF Paul traps are excellently approximated as harmonic. 
Nonlinearity in ion motion was observed when several 
ions are trapped due to their mutual Coulomb repulsion. 
Here nonlinearity couples between the ion-crystal normal 
modes, even at the single quantum level [TH[T3]. Trap 
nonlinearities are important in the context of resonance 
ejection in high resolution mass spectrometry jl4j. How- 
ever, in these experiments ions are typically not laser- 
cooled and furthermore the effect of Coulomb nonlinear- 
ities in the large ion cloud is intertwined with that of 
the trap. Recently, amplification saturation of a single- 
ion "phonon laser" , resulting from optical forces that are 
nonlinear in the ions velocity, was demonstrated jl5j . 

Here, we study the nonlinear mechanical response of 
a single laser-cooled ^^Sr"*" ion, in a linear RF-Paul 
trap. The nonlinearity originates from the higher than 
quadrupolar order terms in the trapping potential. We 
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FIG. 1: Schematic diagram of the experimental set-up and 
the relevant energy levels of ***Sr+ ion. The positively biased 
trap end-caps produce a static trapping potential in the axial 
direction with a small anharmonicity. The ion oscillator is 
driven to the nonlinear regime by a small oscillating voltage 
on one of the trap end-caps. Violet and infra-red laser beams 
provide laser cooling and optical pumping. Scattered violet 
photons are collected by an imaging system and directed ei- 
ther to an EMCCD camera or to a PMT. 



find that the ion steady state response is well described by 
the Duffing model with an additional nonlinear damping 
term |16| . Unlike other realizations of nonlinear mechani- 
cal oscillators, both the linear and the nonlinear damping 
components can be precisely controlled. 

Our trap has the canonical linear four rods and two 
end-caps configuration shown in Fig. [T] The distance of 
the ion to the end-caps and rod-electrodes is 0.65 mm and 
0.27 mm respectively. Here we examine only the motion 
along the axial direction of the trap. In this direction, 
trapping is dominated by the static electric potential 
due to a positive constant voltage on the trap end-caps. 
This potential is well approximated to be harmonic with 
ujo/2tt = 438 KHz. However, as the trap end-caps do not 
satisfy the pure electric quadruple boundary condition, a 
small octupolar contribution to the electric field results a 
positive cubic term in the restoring force and an energy 
level difference of fiwo + HAnin where n is the harmonic 
oscillator quantum number and A„;/27r = 0.8 mHz is the 
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nonlinear dispersion. This nonlinearity becomes increas- 
ingly important with growing oscillation amplitude. The 
ion is driven to the nonlinear regime by adding a small 
oscillating voltage to one of the trap end-caps. The ion is 
Doppler-cooled by scattering photons from a single laser 
beam, slightly red-detuned from the 51/2 ^ Piji tran- 
sition at 422 nm. To prevent population accumulation 
in the 1)3/2 meta-stable level we repump the ion on the 
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Pi/2 transition at 1092 nm. 



We measure the steady-state oscillation amplitude of 
the ion as we slowly scan the drive frequency, w, across 
the harmonic resonance, wq. The scan is from lower to 
higher frequencies (positive sweep) or vice versa (negative 
sweep). Photons scattered during the cooling process 
are collected by an imaging system (N.A. = 0.31), and 
are either directed towards an Electron-Multiplying CCD 
(EMCCD) camera or a Photo Multiplier Tube (PMT). 

We measure the amplitude of motion by taking time- 
averaged images of the ion as shown in Figj2j^a). The 
image is then integrated along the direction perpendicu- 
lar to motion to produce a single curve. Every column in 
Fig. [2]jb) corresponds to a curve produced this way, for 
a positive frequency sweep. As seen, the oscillation am- 
plitude increases as the drive frequency approaches wq, 
continues to increase passed wq, until at a given critical 
drive frequency, Wm, abruptly collapses to a significantly 
lower value. We extract the ion oscillation amplitude 
by fitting the curve to the expected time-averaged po- 
sition distribution. The blue and red lines in Fig. [3] 
are the measured amplitudes for positive and negative 
sweeps respectively. Different curves correspond to dif- 
ferent drive amplitudes. The asymmetry and hysteresis 
as well as the abrupt amplitude changes at specific criti- 
cal drive frequencies are a clear deviation from the driven 
harmonic oscillator response. As expected from a posi- 
tive nonlinearity, the oscillator self-frequency is "pulled" 
to higher values at higher oscillation amplitudes. 

In order to measure the phase difference between the 
ion-oscillator and the driving force, we time-stamp each 
photon measured by the PMT within a single drive pe- 
riod to allocate it with a corresponding drive phase. 
The instantaneous photon scattering rate from the cool- 
ing laser beam is determined by the ions' instantaneous 
velocity through its associated Doppler shift [T^- A 
histogram of the measured photon phases is shown in 
Figj2jc) . A clear sinusoidal oscillation of the photon scat- 
tering rate yields the ion-oscillator phase. The columns 
in Figj2]jd) are photon phase histograms for a positive 
drive frequency sweep. As seen, at the critical frequency, 
Wm, a phase jump of 1.2 radians in the oscillator motion 
accompanies the sudden change in oscillation amplitude. 

Our observations are well accounted for by the Duffing 
oscillator model. The Duffing equation of motion is. 
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FIG. 2; Driven ion-oscillator amplitude and phase, (a) Time- 
averaged ion images taken at various drive frequencies, (b) 
Columns are time-averaged images, integrated along the di- 
rection perpendicular to ion-motion, during a positive fre- 
quency scan, (c) A histogram of the number of photons de- 
tected at different driving force phases, (d) Columns are pho- 
ton phase histograms taken during a positive frequency scan. 
The sohd line is the theoretical phase given by Eq[3] shifted 
by a constant to match the peak in the histograms. 



Here x is the displacement of the ion from its equilib- 
rium position, a is the an- harmonic coefficient, [i is the 
linear damping coefficient and k is the drive amplitude. 
The recoil noise inherent to the spontaneous photon scat- 
tering process, which would appear as a Langevin force 
term, is neglected. An approximate solution to Eqjljcan 
be obtained by the multiple scale method [T]. Here the 
solution has the from x(t) = a(t) cos(a;t — 0) , where a(i) 
is a slowly-varying oscillation amplitude and is the os- 
cillator phase. The steady-state solution for a solves. 
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where cr = a; — cjq is the drive detuning. The steady-state 
solution for a, at a fixed fc, vs. drive frequency is shown 
by the black line in the inset of Figj3] Above a criti- 
cal amplitude Cc, a trifurcates into three solutions. One 
solution with small and one with large amplitude, are sta- 
ble, while the third, intermediate amplitude solution, is 
unstable and is positioned on the state-space separatrix. 
This bistablity persists until the high amplitude solution 
reaches a maximal value, am, at which the drive force is 
overwhelmed by damping and the oscillator is forced into 
a single stable solution. Positive and negative frequency 
scans carry the oscillator into the bistability region along 
different stable attractors leading to the observed hys- 
teresis as illustrated by the arrows in the inset. To com- 
pare with our data we independently measure all the pa- 
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rameters in Eq{T] The driving force amphtude, k, is mea- 
sured by observing ion displacement vs. end-cap voltage, 
ojo is measured via ion response in the linear regime. A 
value of a/Air'^ = 1.24 ± 0.03 • 10^^ Hz^/m^ is measured 
using the observed dependence of a„i on Um — — wq, 
the instability detuning, a„i ~ ^8wof„i,/3a. A value 
of ii/2tt ~ 39.2 ± 0.3 Hz, which result in a quality factor 
Q = 5590, is evaluated using the variation of with the 
drive amplitude, = k/{2fiujQ). The blue and red cir- 
cles in the inset are the measured amplitudes, for positive 
and negative scans respectively, showing good agreement 
with the theoretical curve. 
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FIG. 3: Measured oscillator amplitude vs. drive frequency 
for various driving force amplitudes and both positive (blue) 
and negative (red) scans. The inset shows the Duffing model 
calculation (black solid line) and our measured amplitudes 
(blue circles - positive scan; red circles - negative scan) 

The Duffing oscillator steady state phase is given by, 
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The white solid line in Figj2|^d) shows the theoretical 
phase curve vs. drive frequency for our experimental pa- 
rameters, showing good agreement with our data. 

Linear damping is a very good approximation for most 
mechanical oscillators, as typically dissipation originates 
from coupling of the oscillator to an ohmic bath. Re- 
cently, the contribution of non-linear damping to the mo- 
tion of a nano beam resonator was studied [19] ■ In our 
experiment damping results from the change in radia- 
tion pressure vs. ion velocity. When the laser frequency 
is tuned below the cooling transition, the leading con- 
tribution is indeed linear in the ions' velocity. However, 
as the oscillation amplitude increases or the cooling-laser 
detuning reduced, the effect of damping force terms that 
are nonlinear in the ions' velocity increases [T7]- 

To account for nonlinear damping, Eqjljis modified to 
include a term which is cubic in the oscillator velocity. 



Here 7 is the cubic damping coefficient. The steady-state 
amplitude is now a solution of 119j , 
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When 7 > 0, nonlinear damping acts to effectively in- 
crease dissipation for larger oscillation amplitudes. Un- 
like the linear damping does not increase linearly 
with k but is rather limited by the growing dissipation. 
We find fi and 7 by a maximum likelihood fit of the mea- 
sured Qm vs. k curve to the solution of Eq. [5] It is 
instructive to look at the responsivity, x = 2fj,u!oa/k, in 
order to distinguish linear from nonlinear damping |18j . 
In Fig [4] we plot the measured x for positive scans and 
various drive amplitudes, fc, for two different cooling-laser 
detuning values, S. In Figga) S/2Tr = -420 MHz, 7 = 
and the maximal responsivity is seen to be independent of 
k. In Figgb) (5/27r = -160 MHz, w§7/27r = 0.09±0.002 
/im~^Hz and the maximal responsivity decreases as k in- 
creases. The linear dissipation term, /i, is similar in both 
cases. The solid lines are the solutions of Eqj5] showing 
good agreement with the data. 
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FIG. 4; Calculated and measured responsivity, x ~ 2yLU)oa/k, 
for positive drive scans. Different curves correspond to differ- 
ent drive amplitudes. Due to small drifts in (<100Hz) 
each curve was separately shifted on the frequency axis to fit 
the theoretical curve, (a) Linear damping, the cooling laser 
detuning 5/2% — —420 MHz and 7 = 0. The maximal re- 
sponsivity is independent of drive amplitude k. (b) Nonlinear 
damping, S/2tt = -160 MHz and a;g7/27r = 0.09 fj. m"^Hz. 
The maximal responsivity decreases as k increases. 

We next repeat the measurement of /i and 7 for various 
cooling-laser detunings at a fixed repump-laser frequency 
and lasers intensities. The measured /i and 7 vs. S are 
shown in figures [5ja) and [5jb) respectively. To com- 
pare with the theoretically predicted values we write the 
cooling-laser scattering force, 



Fs{x) = hkcTpp{Sc + kcX, Sr + krx), 
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where fc^/r ^-nd S^/r ^re the wave-vectors and detunings 
of the cooling and repump lasers respectively, F = 27r x 
21 MHz is the spectral linewidth of the level and 
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Pp is the Pi/2 population. The damping coefficients are 
therefore given by the appropriate derivatives, 
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Here m is the ion mass. We calculate pp by numerically 
solving the eight coupled Bloch equations, corresponding 
to the population in all states in the 6*1/2, P1/2 and 
levels coupled by the cooling and repump lasers. The 
cubic damping coefficient is highly sensitive to different 
laser parameters due to the presence of dark resonances. 
The solid lines in figures [sja) andjsjb) are the calculated 
p and 7 showing good agreement with our measured val- 
ues. The two lasers intensities and the repump-laser de- 
tuning were used as fit parameters, yielding values that 
agree within 20% with their measured value. 
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FIG. 5; (a) Linear and (b) cubic damping coefficients for 
various cooling-beam detunings. Filled circles are measured 
values and solid lines are calculated using equations [6] and [7| 

An additional nonlinear damping term, proportional 
to x'^x, results from the laser beam finite size (100 ^m 
FWHM) and has an identical effect to that of 7 on the 
steady-state motion [121 IIH] ■ This term was calculated 
to be small relative to 7 [20] and was taken into account 
in FigjSj 

In conclusion, we have driven a single-ion oscillator to 
the nonlinear regime. The ion steady-state motion, show- 
ing a bifurcation into two stable attractors and hysteresis, 
is well described by the Duffing oscillator model with an 
additional nonlinear damping term. Unlike previously 
studied nonlinear mechanical oscillators, here both the 
linear and nonlinear parts of dissipation can be tuned 
with the cooling laser parameters. 

The study of the nonlinear motion of trapped laser- 
cooled ions opens several exciting research avenues. Since 
trapped atomic-ions can be cooled to the quantum 
ground state, they are an excellent platform to study 
nonlinear behavior in the quantum regime. As shown 
in , unlike the simple harmonic oscillator, a Duffing 
oscillator will demonstrate a clear quantum-to-classical 
transition even when classically driven. Moreover, as the 



ion-spin can be entangled with its motion, it will be pos- 
sible to form a coherent superposition of the two attrac- 
tors states of motion. Laser-cooling of a nonlinear driven 
ion-oscillator has several interesting aspects that can be 
further explored. Since the Doppler shifts associated 
with the oscillation amplitudes in the nonlinear regime 
are significant compared with the cooling transition line- 
width, the laser-cooling force is largely nonlinear in the 
oscillator velocity. Furthermore, the thermal state gen- 
erated by laser-cooling is the result of balance between 
the damping force and the inherent heating due to the 
recoil noise from spontaneous photon scattering. Close 
to the Duffing instability, the ion-oscillator response to 
noise is quadrature dependent. One noise quadrature is 
largely enhanced whereas the other quadrature is sup- 
pressed . Laser-cooling in this case is likely to produce 
squeezed states of motion. 
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